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Abstract. We present a detailed Monte Carlo simulation of the population of the hot 
branch of the white dwarf luminosity function. We used the most up-to-date stellar evo- 
lutionary models and we implemented a full description of the observational selection 
biases. Our theoretical results are compared with the luminosity function of hot white 
dwarfs obtained from the Sloan Digital Sky Survey (SDSS), for both DA and non-DA 
white dwarfs. For non-DA white dwarfs we find an excellent agreement with the obser- 
vational data, while for DA white dwarfs our simulations show some discrepancies with 
the observations for the brightest luminosity bins, those corresponding to L > 10 Lq. 



1. Introduction 



Among other applications — see, for instance, the recent review of lAlthaus etaP(l201Ch 



— the luminosity function and space density of white dwarfs provide interesting con- 
straints on the local star formation rate and history of the Galactic disk in the Solar 
neighborhood. However, to address these questions a large white dwarf sample of 
known completeness is required. The SDSS has provided us with such a sample, and 
from it a r eliable white dwarf lum inosity function for hot white dwarfs has been recently 
obtained dKrzesinski et al.ll2009t) . This white dwarf luminosity function has the inter- 
esting particularity that it has been obtained using solely spectroscopically-confirmed 
white dwarfs from the SDSS DR4, and thus constitutes an excellent testbed to check 
not only the white dwarf cooling sequences at high luminosities, but also our ability 
to model reliably some Galactic inputs necessary to compute the luminosity function. 
Here we describe the results of a comprehensive set of Monte Carlo simulations aimed 
to model the hot part of the white dwarf luminosity function for both DA and non-DA 
stars. 



1 



2 



Torres et al. 




2 0-2 2 -2 

log(L/L Q ) log(L/L ) 

Figure 1 . White dwarf luminosity functions for hot DA and DB white dwarfs for 
two di fferent models. The solid line is the luminosity function of iRrzesinski et all 
(2009), while the dashed line shows the result of our simulations. 



2. The Monte Carlo simulator 



We simulated a synthetic population of disk white dwarfs in the Solar neighborhood in a 
sphere of 3 kpc. An extensive description of our Monte Carlo simulator can be found i n 
previous papers (IGarcfa-Berro et al. 1999; iTorres et al. 2002; Garcfa-Berro et al. 1120041) . 
Thus, here we only summarize the most important inputs. We adopted a dis k age of 
10.5 Gyr, a constant star formation rate and a standard initial mass function (IKroupa 
2001). " Velocities were obtained taking into account the differential rotation of the 
Galaxy, the peculiar velocity of the Sun and a dispersion law which depends on the 
Galactic scale height. A double exponential profile was used, with a scale heigth 
h = 250 pc a nd a scale length L = 1.3 kpc. Also, the initial-final relationship of 
Catalan et all (120081) was adopted. The cooling sequences employed depend on the 
mass of the white dwarf, Mwd- If Wwd ^ 1.1 Mq a CO core was adopted, while 
if Mwd > 1.1 Mq an ONe core was used. In the cas e of CO whit e dwar fs with H- 
rich envelopes we used the evolutionary calculat ions o f Rened o et al.l (120101) . while for 
white dwarfs with ONe cores we used those of Althaus et al. (200 7|). For H-d eficient 
white dwarfs we used the cooling sequences of Benvenuto & Alt hausl (| 19971). which 



corres pond to pure He atmospheres, and the bolometric corrections of Ber geron et al. 



(1995). Finally, we used the same selection criteria employed to cull the observational 
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Figure 2. White dwarf luminosity functions for hot D A white dwarfs and diff erent 
assumptions. The solid line is the luminosity function of lRrzesinski e t al. (2009|), the 
dashed line shows our fiducial model, while the dotted line represents the luminosity 
function when NLTE corrections and metal contamination is considered. 



sample of iKrzesinski et all (l2009h . Specifically, only white dwarfs with g > 14 and 
a fully de-reddened magnitude g Q < 19 were selected, whereas the color cuts were 
-1.5 < (u - g) Q < and -1.5 < (g - r) Q < 0. 



3. Results 

In Fig. 1 we show several model luminosity functions for hot DA and DB white dwarfs 
for two different assumptions. In the canonical model we adopt a DA/DB ratio /da/db = 
0.80, independently of the effective tempera ture. In a second s imulation /da/db de- 
pends on r e ff, as obtained from the SDSS (Krzesinski et al.ll2009l) . The solid lines rep- 
resent the luminosity function of IKrzesinski et all (120091) . while the dashed lines show 
our simulated luminosity functions. The upper panels of Fig. 1 show that the luminosity 
function of hot DAs barely depends on the DA/DB ratio, whilst that of non-DA white 
dwarfs depends sensitively on it. Also, it is clear that the agreement with the observa- 
tional data is excellent for the model in which the DA/DB ratio depends on r e ff. Thus, 
we adopt this model as our fiducial one. 

One puzzling characteristic of the observational luminosity function of hot white 
dwarfs is the existence of a plateau at luminosities around log(L/LQ) « 1. In a first 
attempt to explore its origin, we analyzed the effects of NLTE correc tions and of metal- 
licity. We adopted the NLTE corrections of iNapiwotzki et ail d 1999b for temperatures 
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30, 000 K < T eS < 100, 000 K and gravities 6.50 < logg < 9.75. Additionally we 
assumed that 50% of hot DA white dwarfs have metals in their atmospheres. This im- 
plies that their effective temperature is overestimated by 20% at 80 000 K and 0% at 
40, 000 K. The results are shown in Fig. 2. The dashed line corresponds to our fiducial 
model and the dotted line to the model in which both NLTE corrections and metal con- 
tamination have been included. Clearly, these additional effects have a small impact on 
the computed luminosity function. 

Finally, to assess if the plateau of the luminosity function of hot DAs could cor- 
respond to a recent burst of star formation (INoh & Scaloll 19901) . we computed several 
luminosity functions with different burst strengths occurring at various ages. How- 
ever, we found that for the approximate age of a typical white dwarf of 0.6 Mq at 
log(L/L0) « 1, the burst should have occurred very recently. Thus, there are only two 
possible alternatives left. Either the theoretical cooling models are not entirely reli- 
able at these luminosities — possibly influenced by the initial conditions — or, instead, 
since at this luminosity the number of objects is relatively small the observational error 
bars inherent to the (l/V max ) method have been underestimated ( Geijo et ap2006l) . and 
the plateau corresponds to a statistical fluctuation. 



4. Conclusions 

We presented a set of simulations of the luminosity function of hot white dwarfs. Our 
results are in excellent agreement with the observations for DB stars, while the plateau 
of the luminosity function of DAs at log(L/LQ) « 1 cannot be reproduced. At these 
luminosities NLTE effects and metal contamination play a minor role in shaping the 
luminosity function, and the plateau cannot be explained by a burst of star formation. 
Thus, it could be due to a failure of the cooling models or to a statistical fluctuation. 
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